Do we need to know the temperature in prestellar cores? 
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ABSTRACT 

Molecular line observations of starless (prestellar) cores combined with a 
chemical evolution modeling and radiative transfer calculations are a powerful 
tool to study the earliest stages of star formation. However, conclusions drawn 
from such a modeling may noticeably depend on the assumed thermal structure 
of the cores. The assumption of isothermality, which may work well in chemo- 
dynamical studies, becomes a critical factor in molecular line formation simula- 
tions. We argue that even small temperature variations, which are likely to exist 
in starless cores, can have a non-negligible effect on the interpretation of molec- 
ular line data and derived core properties. In particular, "chemically pristine" 
isothermal cores (low depletion) can have centrally peaked C^^O and C^^S radial 
intensity profiles, while having ring-like intensity distributions in models with 
a colder center and/or warmer envelope assuming the same underlying chemi- 
cal structure. Therefore, derived molecular abundances based on oversimplified 
thermal models may lead to a mis-interpretation of the line data. 

Subject headings: astrochemistry — line: profiles — stars: formation — ISM: 
molecules 
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Introduction 



The cooling in starless (or prestellar) cores is believed to be very effective, so it is cus- 
tomary to model these cores as isothermal at the temperature of around 9-10 K. In general, 
this assumption is confirmed both observationally and theoretically. However, different au- 
thors disagree on possible minor deviations from the uniform temperature, which may exist 
in these objects. 

More precisely, various attempts to model the gas thermal balance in starless cores 
basically reach the same conclusion that the gas an d dust ternperat u res are well coupled 



in the core center and have values of 5-7 K (e. g. iGoldsmith 



2001; Zucconi et al. 2001 



Galli et al.ll2002l : iLesaffre et al.ll2005l ). At a low density periphery, heating processes due to 
cosmic rays and some UV irradiation take over, and the gas temperature rises to 15 K and 
higher (depending on the UV field strength), at the same time decoupling from the dust 
temperature. 

On the other hand, recent obse rvational evid e nce fo r th e temperature struc ture in star- 
less cores is somewhat controversial. iTafalla et al.l (120021 ) and iTaf alia et al.l (120041 ) have found 
an isothermal structure at around 9 K for a samp le of starless cores, including the L1544 and 
L1498 cores. Similarly, Ivan der Tak et al.l (120051 ) have found from H2D"'' observations that 
the gas temperature is nearly constant or at least does no t drop inwards significantly in 
L1544. High-resolution observations have led ICrapsi et al.l (120071 ) to the conclusion that 
the gas temperature does go down in the very center of the L1544 core. A similar thermal 
structure (a few K i n the center and mo re than 10 K in the en velope) is inferred for th e 
L183 starless core by iPagani et al.l (120071 ) and for B68 globule bv lPineda fc BenschI (120071). 



Quite surprisingly, a reversed temperature radial profile for B68 is found by iBergin et al. 



( I2OO6I ). These authors argue that the envelope of this globule is actually colder than its 
central region. 

These deviations from isothermality would have only a subtle effect on results of dy- 
namical and chemical modeling. However, there is an aspect, which is more sensitive to 
temperature variations. In order to interpret observational data in terms of abundances 
and velocities, radiation transfer in molecular lines is commonly utilized. The synthetic line 
profiles depend critically on the radial distribution of the excitation temperature, which is 
a non-linear function of the kinetic temperature T. The importance of adequate t empera- 
ture understanding has already been recognized, e. g., by lDi Francesco et al.l (120021 ). These 
authors found that utilization of the non-uniform temperature profile with decreasing T at 
smaller radii increases abundances of C^^O, CS, and H2CO, inferred for B68, by factors of 
2-3 in comparison with the isothermal case. 
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In this Letter, we study the influence of a non-uniform T distribution on synthetic inte- 
grated intensity Jint radial profiles for various CO and CS (and their isotopomers) transitions, 
using the L1544 and L15 21E starless cores as templ ates. Observational data are taken from 
Tafalla et al.l ((2002) and iTafalla fc Santiago! (120041 ) . The radiation transfer model as well 
as the assumed pre-stellar core chemical and thermal structure are briefly outlined in § 2. 
Results for the detailed chemical modeling are presented in § 3. Influence of the thermal 
structure on the inferred analytical abundance representations is discussed in § 4. Results 
are summarized in § 5. 



2. Model 



Synthetic radial proflles of the integrated intens ities in a prestellar core ar e comp uted 
with the Monte-Carlo code URAN(IA), described by IPavlvuchenkov fc ShustoyI (12004r). We 



utilize den sity distributions for the L1544 and L1521E cores, represented in iTafalla et al. 
(l2002h and lXafalla fc Santiagol boo4 ] as n(r) = no/(l + (r/ro)P), with p = 2.4, tq = 2800 AU 



no 



1.35 X 10^ cm-3 for L1544 and p = 2, tq = 4200 AU, uq 



2.7 X 10^ cm-3 for L1521E. 



A model core is assumed to be static, with a uniform radial microturbulent velocity Vt of 
140 m s~^. We tried other Vt values as well and found that our conclusions do not depend 
on particular assumptions about this quantity. 

Three different temperature proflles are considered (Figured]). In models, designated 
'uni', T is assumed to be 8.75 K in L1544-based runs and 10 K in L1521E-based runs at all 
radii. In 'warm' models T decreases from 15 K at the core surface to about 5 K at r = 
(the core center). The 'cold' temperature proflle is flat with T = 8.75 K (L1544) or 10 K 
(L1521E) at r/routcr > 0.2 and linearly decreases at smaller radii to 4 K at the center. 

Two different approaches for estimating the chemical structure of a starless core are 
employed. In § 3 the CO and CS abundances are obtai ned through a detailed chem ical 
modeling of a core. We use the same chemical model as in iPavlyuchenkov et al.l (120061 ) and 
the density proflle for the L1544 core. The core is assumed to be isothermal (T = 10 K). We 
allow the chemistry to evolve for 2 Myr and then use the resulting molecular abundances 
at various times to model the emergent line proflles. Because of the assumption of uniform 
T = 10 K in the chemical models, they are formally inconsistent with the radiation transfer 
models. However, we calculated the same chemical models with the above non-uniform 
temperature proflles, and confirmed that the chemistry is not affected by these relatively 
small T variations. 



In § 4, we use the density profiles for L1544 and L1521E and approximate the ra- 
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dial abundance profiles with pre-defined analytic expressions. Only C^^O(l-O) (L1544) and 
C^^0(2-l) (L1521E) transitions are considered in that section. 



3. Line intensities with complete chemistry 



In the chemical modeling, we use two parameter sets. In the first set we assume a 
low sticking probability of S* = 0.3 (low depletion case). The model core is illuminated by 
interstellar UV field with G = 0.2 and penetrated by cosmic rays with an ionization rate of 
C = 3 X 10~^^ s~^. This case is int ended to re present "chernically pristine" starless cores like 
CB17 (jPavlyuchenkov et al.ll2006l ) or L1521E (IHirota et al.ll2002l ). The other parameter set is 
more appropriate for "chemically evolved" cores (like L1544 or L1498). In this high depletion 
case, we adopt S = 1, ( = 3 x 10~^^ s~^ and the same UV field. Integrated intensity profiles 
for the considered transitions, temperatures, an d chemical paramet ers, convolved with the 
antenna beams, corresponding to observations in lTafalla et al.l (120021 ). are shown in Figure [2J 
The CO abundance in the low depletion case is not significantly time-dependent, so results 
after 2 Myr of evolution are shown. In the other plots, time steps are chosen, which in the 
'uni' model produce results close to observational data for L1544 (ITafalla et al.ll2002l ). 



In the high depletion case, all molecules are heavily depleted (X(max)/X(r = 0) ~ 200 
for CO), and it is of no surprise that the 'uni' and 'cold' models do not differ from each other. 
Both optically thin (C^^O, C^^S) and optically thick lines are generated at intermediate radii 
and are not affected by the temperature in the center. However, a warmer envelope ('warm' 
model) results in a brighter emission, with the effect being the strongest for the optically 
thick CS (2-1) transition. This may have a direct impact on the analysis of observational 
data. If we use the adopted chemical model to interpret this emission erroneously as arising 
from an isothermal object, we would overestimate the core age by a factor of 1.5. 



In the low depletion case with uniform T (Figure [21 bottom left), the central hole in the 
C^^O integrated intensity distribution does not show up at all. Even though the depletion is 
still quite significant (X(max)/X(r = 0) ~ 60), the absolute molecule number density peaks 
toward the center. Also, the uniform T results in a relatively high excitation temperature. 
However, the central depression in Jint readily appears if we adopt a non-uniform temperature 
profile, with T decreasing toward the center or increasing in the envelope. In the former case, 
the ring-like intensity profile arises because of the weaker central emission, in the latter case 
it is caused by the brighter envelope. 

The temperature dependence of the C^^S and CS profiles is much weaker, because these 
profiles are mainly formed in the core envelope. As in the high depletion case, the isothermal 
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model and the model with a cold center do not differ from each other, while the emission in 
the warmer envelope model is brighter than in other models. 

Thus, whether or not the moderate depletion is visible as a central hole in the C^'^O 
(1-0) Jint distribution, depends on the adopted thermal structure. A chemically evolved 
core always looks like a chemically evolved core, while a less evolved core may appear to be 
chemically mature just because its temperature structure is not taken into account properly. 

One may suggest that the temperature structure can be inferred from transition ratios, 
but in this case the interpretation is also not straightforward. In Figure [3] (top panel) we 
show the T profiles which would have been derived from the ratios of integrated intensities 
C^^O(2-1)/C^^O(1-0), assuming these transitions are thermalized and optically thin. Cores 
with warmer envelopes do stand out on this plot, but the C^^O(2-1)/(1-0) ratios do not 
allow to distinguish between 'cold' and 'uni' models. In particular, the central temperature 
drop in 'cold' models is only formally present. 

Some hints can, possibly, be extracted from the analysis of higher transitions. On the 
bottom panel of Figure [3] we show the computed C^'^0(3-2)/C^^O(2-l) integrated intensity 
ratios. In the pristine models the (2-1) line is saturated (i.e., its optical depth is greater 
than 1), while the (3-2) transition is (moderately) optically thin. Because of that, higher 
density at the core center enhances the (3-2) transition relative to the (2-1) transition. 
Thus, the C^^O(3-2)/C^^O(2-l) ratio is higher at small r in 'uni' and 'cold' pristine models. 
In 'warm' models the computed 0^^0(3-2)70^*^0(2-1) ratio is large because in the warm 
envelope collisions more effectively populate higher levels. 

In 'uni' and 'cold' evolved models , both 0^^0(2-1) and 0^^0(3-2) lines are optically 
thin, and the (2-1) transition does not show the saturation effects. At the same time, the 
(3-2) line is weak due to low density at the periphery and significant depletion at the core 
center. Thus, the 0^^0(3-2)/0^^0(2-l) ratio is low. 

Note the significant difference in 0^^0(3-2)/0^^0(2-l) ratios between the 'cold' and 
'uni' pristine models in the low depletion case (Figure [31 bottom panel, solid and long- 
dashed curves). This implies that in surveys of "chemically young" starless cores, which are 
currently underway, it is still desirable to pay specific attention to higher transitions in order 
to assess the thermal structure of these objects. 

Another issue, which has been brought to our attention by the referee, is related to 
isotope selective photodissociation, which preferentially destroys O^^O. This process is not 
taken into account in our model, but it may affect the ability of O^^O to trace warm envelopes 
of starless cores. From this point of view, ^■^OO transitions may be better envelope probes. 
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4. Analytical fits 

To study the sensitivity of derived molecular abundances on the ter nperature structure , 



we applied t he technique, utilized by Tafalla and co-authors for L1544 (ITafalla et al.ll2002l ) 



and L1521E (ITafalla fc Santiago! l2004l ). i. e., we approximate the molecular abundance pro- 
files with X{r) = X^exp [— n(r)/nd] and combine them with the three T profiles described 
above. 

LI544 — L1544 is a well known example of a starless core with some species, like CO, CS, 
C2S, being strongly depleted. Thus, according to the results of § 3, we do not expect details 
of the thermal structure to have any significant effect on the interpretation of observational 
data. 

With each of the three temperature profiles, we computed a series of RT models for 
different values of Xq and rid and compared the resultant distributions of the integrated 
intensity with the observed distributions. In FigurelHwe present the results of the comparison 
as a function of Xq and Ud- The best-fit positions in each plot are indicated by triangles. 
Combinations that bracket the scatter in observational points for the C^^O(l-O) transition 
are shown with crosses and boxes. The error bars in Figure H] represent the ext ent of scatter 



in th e observed integrated emission as a function of core radius (e.g. Fig. 3 in ITafalla et al. 



2OO2I). 



All the considered models do formally agree on the inferred values of Xq ^ 10~^ and 
TT-d = 8 X 10^ cm~^. However, the overall range of these parameters allowed by the scatter 
in Jint differs for various models. It appears that observations are marginally consistent 
with abundance profiles in which CO molecules are either totally depleted (crosses) or very 
moderately depleted (boxes). As the scatter mostly reflects the complex internal structure of 
the core, a difference between the bracketing Xq — Ud combinations may, to a certain degree, 
serve as a quantitative manifestation of a chemical inhomogeneity of the core. 

L1521E — Unlike L1544, the starless core L1521 E seems to be rich e r in gas phase CO, 



which is interpreted as a sign of its chemical youth (IHirota et al.l l2002l : ITafalla fc Santiago 



20041 ). Indeed, if we use the same exponential X{r) profile as before and a uniform T of 10 K, 
we find the best agreement with observations of the 0^*^0(2-1) transition for Xq = 10^'' and 
Ud = 10® cm~^. This value is greater than the central density uq of L1521E, i. e. the 
abundance distribution is near ly uniform. This is why in the 'uni' model our conclusions are 



identical to those obtained by ITafalla fc Santiago! (120041 ). 



However, the fit is different if we assume 'cold' or 'warm' T profiles. Neither centrally 
depleted nor constant abundances are able to reproduce the observed Jint in L1521E. In 
order to check other possible inferences for the L1521E chemical composition, we employed 



- 7- 



a different X{r) profile as well, assuming that X{r) = Xq exp(— r/rco)- In Figure [5] we com- 
pare the fittin g results for the C^^0(2-l ) transition in the L1521E core to the observational 
data, given in iTafalla &: Santiago! (120041 ) . A combination of uniform T and nearly constant 
abundance is able to fit the observations almost perfectly (solid line). To reproduce the 
observed Jint{f) in a core with a non- uniform T distribution, one need a centrally peaked 
abundance profile (Figure [5]). Again, like in a non-evolved chemical model, in models with 
lower central T we may get either flat or even depressed Jint('") despite the fact that the un- 
derlying abundance is growing toward the core center. The fit quality is not high, especially, 
in the 'warm' model, but this problem can be solved by using an abundance profile with 
central flattening (like exp[— r^/rg]) or by constructing an adequate chemical model. 



5. Conclusions 

The idea of this study is to draw attention to the problems related to the thermal 
structure of prestellar cores. We demonstrate that realistic radial temperature profiles can 
produce quite different line intensities, even when the underlying chemical structure is the 
same. Effects of a non-uniform T distribution can mimic effects of chemical evolution and, 
thus, may hamper the interpretation of observations. This is especially true for "pristine" 
starless cores and is less important for chemically evolved cores. Finally, we would like 
to note that in this study we only consider transitions of CO and CS molecules. While 
ammonia observations represent a more conventional tool to measure the temperature in the 
interstellar medium, they are, probably, prone to the same kind of uncertainty. 

This study is supported by the DFG grant "Molecular Cloud Chemistry" HE 1935/21- 
1. DW acknowledges support from the RFFI Grant 07-02-01031 and the Russian Science 
Support Foundation grant. 
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Fig. 1. — Model temperature profiles. In the 'uni' profile case, the temperature is uniform 
across the core at 8.75 K for L1544 and 10 K for L1521E (not shown). The 'cold' profile 
corresponds to the cold center, while in the 'warm' case, the envelope is warmer than in the 
'uni' case. 
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Fig. 2. — Integrated intensity profiles for C^^O(l-O), C^^S(2-1), and CS(2-1) transitions, 
calculated with different assumptions on the core thermal structure. Thick grey lines and 
right axes represent the relative abundances of C^^O and CS. Shown are results for the high 
depletion case (top row) and the low depletion case (bottom row). Time moments are chosen, 
which provide the best agreement with L1544 observations in the 'uni' model, except for the 
bottom left panel, where results at the end of the computation are used (see text). 
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Fig. 3. — Top: Temperatures computed from C^^O(2-1)/C^^O(1-0) integrated intensity 
ratios, assuming these transitions are thermalized and optically thin. Bottom: 0^*^0(3- 
2)/C^^O(2-l) intensity ratios. Abundances for t = 2 Myrs are used. 
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Fig. 4. — Results of the analytical fits for 'uni' (left column), 'cold' (middle column), and 
'warm' (right column) models. Top: The observed (connected solid circles) and computed Jmt 
profiles for the best-fit Xo — rid combination (triangles) as well as for bracketing combinations 
marked with crosses and boxes on bottom panels. Bottom: Xq — diagrams indicating the 
quality of agreement with observations. The darkest color corresponds to the best agreement. 
The best-fit positions for each T profile are also marked with triangles. 
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Fig. 5. — Results of the analytical fits for the L1521E core. Radial abundance profiles which 
produce the best fit to the data are shown in the right panel with solid line ('uni' model), 
dashed line ('warm' model), and dotted line ('cold' model). Corresponding J^^t profiles 
are shown in the left pan el with the same line styles. Observational data are taken from 



Tafalla fc Santiagol (12004] ). 



